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1. Introduction 

Investigation of the geometry of the set of quantum states remains a subject of current scientific 
interests in view of possible apphcations in the theory of quantum information processing. 
The set Qjy of pure quantum states belonging to a dimensional complex Hilbert space 
7^ AT is known to be equivalent to the complex projective space, Ojv = <CP^~^, of — 2 
real dimensions. However, as this set is embedded into the A^^ — 1 dimensional set Qat of 
density matrices of size by a non- linear constraint, p = p"^, the geometric structure of the 
set of mixed quantum states is rather involved [HE]- The only simple case corresponds to the 
one-qubit system, N = 2. 

The set ^2 of N = 2 pure states forms the Bloch sphere, <CP^ = S'^, with respect to the 
standard Hilbert-Schmidt metric. The 3-disk inside the sphere, often called the Bloch ball, 
represents the set Q2 of one-qubit mixed states. In this simple case any projection of this set 
onto a plane forms an ellipse, which can be degenerated to an interval. In the case of A'^ = 3 the 
8 dimensional set Q3 of one-qutrit mixed states is neither a polytope nor an ellipsoid [3llH[5], 
and the set ^3 = CP^ of its extremal states is connected and has four real dimensions. 

Due to the high dimensionality of the problem our understanding of the geometry of the 
set Qat of mixed states is still rather limited. This set forms a convex body which contains 
an in-ball of radius r^v = \Jl/N{N — 1) and can be inscribed into an out-sphere of radius 
Rn = {N — l)rp4 = ^ {N — l)/N [2j. Some information on the subject can be gained by 
studying the 2-dimensional cross-sections of Qn as demonstrated in [6l [71 [8] for A^ = 3 and 
A^ = 4. Another option is to investigate projections of this set into a plane - such an approach 
was advocated for A^ = 3 in [9]. As the set Qat of quantum states is convex, also its cross- 
sections and projections inherit convexity. 

In this work we study the general structure of a two-dimensional projection of the set 
Qn of mixed states. A bridge between the geometry of the set of quantum states and the 
notion of numerical range used in operator theory is established. For any operator A, acting 
on the complex Hilbert space TIn, one defines its numerical range \10\ [TT] (also called field 
of values) as a subset of the complex plane which contains expectation values of A among 
arbitrary normalized pure states, 

W{A) = {z:z = {i^\A\i;), \,p) G Hn, (V'I^) = !}• (1) 

We analyse the set of orthogonal projections of the set Qn onto a 2-plane and prove 
that it is equivalent to the set of all possible numerical ranges of complex matrices of order A^. 
Numerical ranges of normal matrices of size A^ correspond to orthogonal projections of the set 
Cat of classical states - the (A^ — l)-dimensional simplex A7V-1 C M^""*^. 

Further information on the structure of the set of quantum states of a size A^ can be 
obtained by studying the numerical shadow |12l \TE[ [Ti] of various matrices of order A^. For 
any operator A acting on one defines a probability distribution Pa{z) on the complex 
plane, supported in the numerical range Ty(A), 

Pa{z):= [ dp{^)5(z-{i;\A\i;)). (2) 

Here p{tp) denotes the unique unitarily invariant (Fubini-Study) measure on the set Oat of 
A'^-dimensional pure quantum states. In other words the shadow P of matrix ^4 at a given 
point z characterizes the likelihood that the expectation value of A among a random pure 
state is equal to z. 

The distribution Pa{z) is naturally associated a given matrix A, and some of its properties 
were described in [13j. In this work we advocate a complementary approach and show that 
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investigating the shadows of several different complex matrices ^4 of a fixed size N contributes 
to our understanding of the structure of the entire set Qat of quantum states. In a sense, the 
choice of a matrix A corresponds to the selection of the plane, onto which the set of quantum 
states is projected. 

This paper is organized as follows. In section[2]we fix the notation and introduce necessary 
concepts. A link between 2-dimensional projections of the set of quantum states of a given size 
N and the set of possible numerical ranges of matrices of order is presented in section |3j 
In section |4] we analyze different classes of numerical shadows of matrices of small order 
N = 2, 3, 4. Unitary dynamics of a pure quantum state in the background of numerical shadow 
is presented in section [5} Section [6] is devoted to the mixed-states numerical shadow, which 
corresponds to a projection of the full set Qat of density matrices onto a plane. The case of a 
large dimension, ^ 1, is treated in section [7] jointly with the shadow of random matrices. 
Finally, in section |8] we provide some concluding remarks and summarize the contribution of 
this paper. 

2. Classical and quantum states 

Let p = {xi,X2, . . . , xat} be a normalized probability vector, so j;^ > and J2iLi — 1- Such 
a vector represents a classical state, and the set Cat of all classical states forms an (A — 1) 
dimensional regular simplex Ajv-i C R^""*^. There exist exactly A^ classical pure states, which 
correspond to the corners of the simplex. All other classical states can be expressed by a convex 
combination of pure states and are called mixed. Typical mixed states are characterized by 
the full rank and they form the entire interior of the probability simplex. 

In quantum theory one describes a physical system with A distinguishable states by 
elements of a complex Hilbert space T-Ln of size N. Its elements represent pure quantum 
states, 1-0) G T~tN- Quantum states are assumed to be normalized, HV'lP = (V'lV') = 1; so they 
belong to the sphere of dimension 2 A — 1. Since one identifies two states, which differ by a 
global phase only, |^) ~ = e~°'\ip), the set of all pure quantum states Qjy, which act on 
7^ AT, is equivalent to the complex projective space, 0,]\f = (CP'^~^ p]. 

In analogy to the classical case, one also defines mixed quantum states (density matrices) 
by a convex combination of projectors onto pure states, p = where Pi > and 

J2iPi — 1- Let us denote the set of all density matrices of order A^ by Qa^- It contains all 
density operators which are positive and normalized, 

Qn = {p:'Hn^ nN, P* = P, P>0, TTp = l} . (3) 

Since density operators are Hermitian and normalized, this set is A^ — 1 dimensional. It 
includes the set of classical states, Qn D Cn = Aat-i, as well as the set of pure quantum 
states, Qj^Dnj^ = GP^'K We are going to work with the geometry implied by the Hilbert- 
Schmidt norm of a matrix, |^|hs ■= y^Tr(^*A), and the Hilbert-Schmidt distance in the space 
of matrices, 

dns{A, B) := \A - B\ns = VT^{A -B){A- B)*. (4) 
It will be also convenient to define a real inner-product by setting the polar identity 

{A, B) = \\A + B\ls-\\A-B\l^ = l[Tr{A*B + B*A)]. (5) 

If ^* = ^ and B* = B then {A, B) = ivAB. 

In the set Oa^ of quantum pure states one defines the Fubini-Study measure //pS) which is 
induced by the Haar measure on U {N) and is invariant with respect to unitary transformations. 
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In the case of one-qubit states this measure corresponds to the uniform distribution of points 
on the Bloch sphere S'^. 

In practice, to generate pure states at random according to the measure fips it is sufficient 
to generate uniformly points at the sphere S*^-^"^. One may also select an arbitrary column, (or 
row) of a random unitary matrix U distributed according to the Haar measure. It directly gives 
the set of N coefficients of the random state in a given basis, = X^^iCj|i). For instance, 
choosing the first column of U we set Cj = Ui^i for i = 1,...,N. Alternatively, one may 
generate N independent complex random numbers Zi and renormalise them, a = Zij \/X]t 
to obtain the desired distribution [151 E] • 

In this work we are going to use the following 

Proposition 1 Let \'\\)) E a random pure state of size N distributed according to the 

Fubini-Study measure. If one represents it in an arbitrary fixed basis, = YliLiCi\i) then 
the squared absolute values of the coefficients, pi = |cip, form a probability vector distributed 
uniformly in the probability simplex Aat.i. 

This is equivalent to the known statement (see e.g. [2|), that the only constraint on the 
components of a single column of a random unitary matrix U distributed according to the 
Haar measure is the normalization condition, P{Uu, . . . Uni) ~ (5^1 — |C/jiP^- This fact 
directly implies 

Corollary 2 For any quantum state p define a classical state p = diag(/9), so pi = pa. Then 
the Fubini-Study measure on the set Vt]\j of quantum pure states induces by this mapping the 
uniform measure in the classical probability simplex Ajv-i- 

In the case of A'^ = 2 the Fubini-Study measure covers uniformly the Bloch sphere S*^. 
Working with the standard polar coordinates, (r, 9, ip), we write the element of the volume of 
the unit sphere as dS = d(ps'm9d9 = d(pd{cos9). The polar angle 9 is defined with respect to 
the axis z, so the projection of a point of the sphere at this axis reads z = cos 9. Hence the 
Fubini-Study measure implies the uniform distribution d{cos 9) = dz along the one-dimensional 
set Ai of = 2 classical states. 

3. Numerical range as a projection of the set of quantum states 

The set = CP^^^ of pure states of size A^ forms the set of extremal points in Q^r. Any 
mixed state p £ Qn can be thus decomposed into a convex mixture of projectors The 
expectation value of an operator A among a pure state reads = TipA. Taking into 

account the convexity of VK(A), the standard definition ([T]) of the numerical range of A can 
be therefore rewritten as [16] 

W{A) = {z:z = TipA, p E Qn}. (6) 

This expression suggests a possible link between numerical range and the structure of the 
set the Qn- Usually one studies numerical range ly(^) for a given A [11]. Here we propose 
to fix the dimension N and consider the set of all possible numerical ranges of matrices A of 
this size to analyze the geometry of quantum states. More precisely, we establish the following 
facts. 

Proposition 3 Let Cn denote the set of classical states of size N , which forms the regular 
simplex A^v-i in 1R,^~^. Then the set of similar images of orthogonal projections of Cn on a 
2-plane is equivalent to the set of all possible numerical ranges W{A) of all normal matrices 
A (such that AA* = A* A) of order N. 
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Proposition 4 Let Qn denote the set of quantum states size N embedded in IR-^^"^ with 
respect to Euclidean geometry induced by Hilbert- Schmidt distance. Then the set of similar 
images of orthogonal projections Qm on a 2-plane is equivalent to the set of all possible 
numerical ranges W{A) of (arbitrary) matrices A of order N . 

To prove the above propositions we will need an abstract lemma concerning the real 
inner-product Euclidean spaces. 

Lemma 5 Suppose ui,U2,vo € V, where V is a Euclidean vector space (with inner product 

1 /2 

(•,•) and norm \x\ = {x,x) ' ), vq and dim (span {ui, U2, vq}) > 2. Then there exist real 
numbers a > 0,71,72 such that the vectors 

vi := - (ni + 71 Wo) , ^2 := - (^2 + J2V0) (7) 
a a 

are normalized and orthogonal, 

= 1 = |f2|^ , {V1,V2)=0. (8) 

Proof. Let u[ = Ui — ^^^^^WO) « = 1,2. By hypothesis + > 0. For i = 1,2 set 
Ci '■= |fo|7i + ^^\vgf so that Vi = ^ (^u'^ + |:^^^oy The desired equations become 

|/|2|2 2 |/|2|2 2 / / /\i n /'n\ 

1^1 1 + Ci = a , \U2\ + C2 = a , {ui,U2) + C1C2 = U. (9j 
Eliminating coefficient a we arrive at a quadratic equation for cf or c^. Set 



then 



Li "'21 ' 



C? =\{\u'2?-\<?) + \^d, (10) 

Cl =\{\u[?-W2?) + \^d, (11) 

sign (C1C2) = -sign (u'l, Ug) , (12) 



a 



1 /I / i9 1 ( lON 1 



+ + (13) 



Recall l^i'iP + |'U2p > by hypothesis thus q > 0. There are generally two solutions 
differing only in the signs of ci and 02- If {u[,U2) = then y/d = — and one of 

the three following cases apply: 

1. \u[\ > Itigl > 0, Cl = 0, C2 = ±Y^mp — a = \u'i\; 

2. IU2I > Iti'il > 0, Cl = zby^mP — C2 = 0, a = IM2I; 

3. \u'i\ = \u2\ > 0, Cl = 0, C2 = 0, a = \u'i\. □ 

Note that the formulae (10, 11) for ci and C2 allow us to obtain the constants 71 and 72, 
which enter eq. Iu\. The scaling factor a = 1 if and only if {u\,U2)'^ = (1 — mp)(l — mP), 
l^'iP < 1 and 1^2^ < 1. 

This lemma implies the following 

Corollary 6 Suppose E C {x £ V : {x, vq) = 1} and ui,U2 £ V define a linear map ^ : E ^ 
C by X ^ {x,ui) + i(x,M2)- Unless ui,U2 G Mwq in which case $ is constant, the map <I> is 
isometrically isomorphic to an orthogonal projection followed by a similarity transformation 
(dilation and translation). 
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Proof. By Lemma [5| there exist orthonormal vectors Vi = -{ui+^iVo) for i = 1,2 and 
a > 0. Let Vq = span {t;i, t;2}. The orthogonal projection onto Vq is given by irx := 
{x,vi)vi + {x,V2)v2 and this is the general form of a rank 2 orthogonal projection. The 
linear map 6 : aivi + 02^2 1— )■ ai + ia2 is an isometry Vq — >• C. If x G E then {vq, x) = 1 and 

2 2 
{oTTx - (71 vi +72^2)) = 0'^{{ui+'yiVo,x)-'yi)vi = 9^^{ui,x) Vi = (/)x.(14) 

i=l i=l 

□ 

Now we are ready to prove the main result of this paper, namely Proposition [3] and 
Proposition |4j 



3.1. Normal matrices 

Proof of Proposition Let ^ be a normal matrix of order with eigenvalues 

{Ai, . . . , Aat}. With respect to an orthonormal basis of eigenvectors of A one has 
Yli^j=i'^i^ij'^j — Si^i'^jlV'iP a-nd the numerical range Wa is the image of the simplex 
AAr_i := |t € : > Vi, J2^^i *i = l} under the map 

N N 

— > ^tjReAi + i^tiImAi = (t, ui) + i(t, 112), (15) 

1=1 i=l 

where t G A7V-1, and (ui)- = ReAj,(n2)j = ImAj for 1 < i < A^. If A ^ cl (multiple 
of the identity, the eigenvalues are all equal) then Lemma [5] and Corollary [6] apply with 
Vq = (1, . . . , 1), which completes the proof of Proposition [sj □ 



3.2. Non-normal matrices 

Proof of Proposition The set Qn of quantum states ([s]) contains Hermitian operators p 
which can be diagonalized, p = UDU*. Here U is unitary while D is a diagonal matrix with 
da > and J2iLi da = I. 

Consider any matrix A of order and write Tr pA = TipAi + iTr/9^2 with Ai = 
^ {A + A*) and A2 = {A — A*). Lemma 5 and Corollary 6 apply now to the map 

$ : p I — > TvpAi + i TvpA2 (16) 

of the set Q^r onto numerical range 1^(A) with V representing the linear space of complex 
matrices of size A'^ (or the real subspace of Hermitian matrices), the real inner product ([5]), 
and Vq = I, ui = Ai, U2 = A2 provided A ^ cl. 

Thus we have shown that for any matrix A its numerical range VF(^) is equal to an 
orthogonal projection of the set of density matrices. To show the converse we may read 
formulae ([7| backwards: the projection of Qat is determined by two orthonormal Hermitian 
matrices Vi and V2, which then satisfy |Vi|hs = |^2|hs = 1 and Tr(FiV2) = 0. Set A = Vi+iV2, 
which gives now the required matrix such that VF(74) is equal to the desired projection. In 
this way a link between numerical ranges of generic matrices of order A^ and projections of 
the set Qn onto a two-plane is established and Proposition [4] is proved. □ 

To obtain explicit formulae for the similarity transformation corresponding to an arbitrary 
matrix A of order A^ define three traceless matrices 

Tr A ^ 1 

B = A-—I, B, = -{B + B*), B2 = -{B-B*). (17) 
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The latter two represent vectors in the Hilbert-Schmidt space and correspond to u\ , U2 in 
Lemma 5. Making use of the Hilbert-Schmidt norm we compute the required coefficients for 
a given traceless matrix B 



4. Numerical shadow and quantum states 

The projectors onto pure states form extremal points of the set Qn of quantum states, 

hence the shape of a projection of the set I^at of pure states onto a given plane coincides with 
the shape of the projection of the set of density matrices on the same plane. As shown in the 
previous section this set is equal to the numerical range of a matrix A of size N, which 

determines the projection. 

However, the differences appear if one studies not only the support of the projection but 
also the corresponding probability measure. A measure -Pa (2) determined by the numerical 
shadow ([2]) is induced by the Fubini-Study measure on the set of pure state. Thus the 
standard numerical shadows of various matrices of size can be interpreted as a projection 
of the complex projective space, JIat = CP^^^, onto a plane. Before discussing in detail the 
cases of low dimensions, let us present here some basic properties of the numerical shadow |13| 
(also called the numerical measure |14j). 

1. By construction the distribution PAiz) is supported on the numerical range of Ty(^) and 
it is normalized, PA{z)d'^z = 1. 

2. The (numerical) shadow is unitarily invariant, Pa{z) = PuAU*{z)- This is a consequence 
of the fact that the integration measure d^{ip) is unitarily invariant. 

3. For any normal operator A acting on Hn, such that AA* = A* A, its shadow covers the 
numerical range Ty(j4) with the probability corresponding to a projection of a regular 
A'^-simplex of classical states Cn (embedded in R^~^) onto a plane. 

4. For a non-normal operator A acting on Hn, its shadow covers the numerical range H^(A) 
with the probability corresponding to an orthogonal projection of the complex projective 
manifold flN = <DP^~^ onto a plane. 

5. For any two operators A and B acting on Ti^, the shadow of their tensor product does 
not depend on the order, 



To show this property define a unitary swap operator S which acts on a composite Hilbert 
space and interchanges the order in the tensor product, S{\z) (8) \y)) = \y) (g) \z). Thus 
{x\A(i^ B\x) = {x\S* B AS\x), and since S is unitary it does not influence the numerical 
shadow induced by the unitarily invariant Fubini-Study measure on complex projective 
space. 



and sign (C1C2) 




PA(g,B{z) = Pb®a{z) ■ 



(20) 
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4.1. One-qubit states, N = 2 

Analysis of the numerical shadow is particularly simple in the case of matrices of order N = 2. 
The spectrum of the operator A consists of two complex numbers, (7{A) = {Ai, A2}. 

In the case of a normal matrix A the numerical range forms the closed interval 

[Ai, A2], and the numerical shadow Pa{z) covers this interval uniformly [13j . 

If the matrix A is non-normal the numerical range forms an elliptical disk with Ai, A2 as 
focal points and minor axis, d = ^^Ti AA* — |Aip — |A2P- For a simple proof of this 1932 result 
of Murnaghan fT7] see the note by Li [18] . In this generic case the numerical shadow is given 
by the probability distribution obtained by the projection of the hollow Bloch sphere of one- 
qubit pure states onto a plane [13]. In particular, the cross-section of the numerical shadow 
supported in an interval x G [0, 1] is given by the arcsine distribution, P{x) = l/{TTy^x{l — x)). 
Non-normal case is shown in Fig{T| obtained for a matrix. 



A 



(2) 



ao 



-1 



For simplicity we have selected the centred matrix such that TrA = so that one has B = A 
in eq. (17). The normalization constant oq = -y/2/5 is chosen in such a way that the scaling 
constant defining the projection in (18) is set to unity, q = 1, so the shadow of the set 



of quantum states is shown in its 'natural size': The distance between both eigenvalues, 
/ = 2ai = -v/2j is equal to the diameter of the Bloch ball, 2i?2 = V^. 



a 0.0 




-0.4 



-0.4 0.0 0.4 
3? 




f 2") 

(a) Shadow of matrix ' (b) Cross-section for 3 = 



Figure 1. Projection of the set ^2 of one-qubit states generated by the numerical 
shadows of operators of order = 2 a) numerical shadow of generic matrix AJ, with 
an elliptical support. Eigenvalues are denoted with crosses and dashed circle of radius 
i?2 — \f^l'^ denotes the diameter of the Bloch ball. Numerically obtained histogram is 
plotted in black, analytical plot is blue. Plot is done for matrix translated in such a way 
that its trace (★) is equal to zero and suitably rescaled. b) Histogram of cross-section 

, solid line represents a probability 

-1 



1 1 

\/2' \/2 



of the shadow supported in the interval 
density function of the arcsine distribution P[x) 



4-2. One-qutrit states, N = 3 

The structure of the numerical range for = 3 was analyzed in detail by Keeler et al. |19j . 
Numerical range of a matrix A of order = 3 with spectrum Ai, A2, A3 forms: 
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a) a compact set of an 'ovular' shape with three eigenvalues in its interior; 

b) a compact set with one flat part - e.g. the convex hull of a cardioid; 

c) a compact set with two flat parts - e.g. the convex hull of an ellipse and a point outside 
it; 

d) triangle. For any normal matrix A its numerical range is equal to the triangle spanned 
be the spectrum, 1^(74) = A(Ai, A2, A3). In the latter case the numerical shadow can be 
verbally interpreted as the shadow of the set C3 of iV = 3 classical states - a uniformly 
covered equilateral triangle A2. 
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Figure 2. Projections of the set of one-qutrit states generated by the numerical 
shadows of operators of order = 3; a) a generic matrix with an oval-like 
numerical shadow, b) J^^^ with one flat part of the boundary dW of the numerical 
range, c) a simple sum with two flat parts of dW , d) a diagonal normal matrix 
Ag^^ with numerical range equal to the triangle of eigenvalues, represented by (+). 
Dashed circle of radius i?3 represents projection of the sphere in which fig is inscribed. 
All plots are done for matrices translated in such a way that their trace (7k-) is equal to 
zero and suitably rescaled. 



The four classes of = 3 numerical ranges are illustrated in figure [2} It shows the 
numerical shadow supported on the corresponding numerical range, obtained for 
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The symbol ojk denotes the k-th root of unity, so ws = exp(i27r/3). As before the matrices 
are chosen to be traceless, so i? = ^ in (17) and the shadows are centered. Furthermore, 
the normahzation constants are designed to assure that the scahng constant in eq. (18) in 
every case is set to unity, a = 1, so the Figure shows images of the set of quantum states 

(3) 

in its natural size. For instance in the case of the diagonal matrix ^3 the prefactor reads 
03 = y^2/3, so that the eigenvalues are located at the distance \/2/3 from the origin. This is 
just the radius R3 of the sphere in which the set Q3 is inscribed. 

The study of the geometry of the numerical range was initiated by Kippenhahn |20| and 
later developed by Fiedler [21] . In more recent papers \22\ [23] the differential topology and 
projection aspects of numerical range were investigated. In particular it was shown [22j that 
the numerical range of a generic matrix A of order three pertains to the class a) above, as 
the boundary of VF(^) does not contain intervals. Critical lines inside the range, analysed in 
|22^ [23] where shown to influence the structure of the numerical shadow [13] . Thus we may 
now relate the critical lines with the geometry of complex projective spaces projected onto a 
plane. 

In the one-qutrit case = 3 obtained probability distributions can be interpreted as 
images of the set of pure states ^3 = CP^ on the plane. Although it is not so simple to 
imagine the structure of the complex projective space [24], some experience is gained by 
studying numerical shadows of various non-normal matrices of size 3. 



4-3. Four-level systems, N = 4 

Various shapes of the numerical range for matrices of size A^ = 4 correspond to various 
projections of the set Q4 of quantum states of size four. As in the case of the qutrit we analyse 
numerical shadows of traceless matrices normalized such that the scaling constant a is set to 
unity. 

Even though several results on geometry of numerical range for A^ = 4 are available 
|25l [26] , complete classification of numerical ranges in this case is still missing. To provide an 
overview of the possible structure of the numerical shadow we analysed the following matrices 
of order four. 
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3i 9i 3i 



(g) (h) (i) 

Figure 3. Projections of the set 5^4 of iV = 4 quantum states emerging as numerical 
shadows of appropriately normalized operators of size 4: a) a generic matrix Aq^"* with 
a an oval-like numerical range W{A), b) A[^^ with one flat part of the boundary dW 
of the numerical range, c) A'^^ being a simple sum 3 1 with two flat parts of dW, 
d) ^3**^ - a simple sum 2 © 2 with two flat parts of dW, e) a''^^ three flat parts of dW 
connected with corners and one oval-like part, f) Ag'^'' three flat parts of dW with only 
one corner and two oval-like parts, g) A'^^ a simple sum 2 ® 2, with four flat parts 
of dW, h) Aj"^^ pair of flat parts of dW connected with a corner connected with two 
oval-like parts, i) a diagonal normal matrix A^'^ with numerical range W equal to the 
convex hull of eigenvalues denoted by (+). All plots are done for matrices translated 
in such a way that their trace (*) is equal to zero and suitably rescaled. 
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Numerical shadows of these representative of each class of = 4 matrices are shown in 
figure[3j The pictures can be interpreted as projections of the 6-dimensional complex projective 
space CP^ onto a plane. Making use of formula (18) we find that the normalization constant 
for the last example ^g^^ reads as = l/\/2. Thus the diameter of the shadow, 2a4 = 
coincides in this case with the Hilbert-Schmidt distance between any two orthogonal pure 
states in Q4. The dashed circle of radius i?4 = \/3/2 represents projection of the sphere into 
which the set Q4 can be inscribed. 



5. Unitary dynamics projected inside the numerical shadow 

As the numerical range and the numerical shadow give us an opportunity to observe the 
structure of the space of quantum states, it is possible to apply these tools to investigate 
quantum dynamics. A unitary time evolution of a quantum system is governed by the 
Hamiltionian operator H (i.e. a self-adjoint operator representing the total energy of the 
system), which leads to U{t) = exp(— ii/t). 

Note that eigenvalues of the Hamiltonian determine the cyclicity of the trajectory. The 
trajectory is periodic iff eigenvalues of the Hamiltonian are commensurable. The period in 
this case is given by least common multiple of the eigenvalues. 

Let us consider a three-level system (qutrit). For concreteness, we choose the 
Hamiltionian 

-1-i 1 



H 



-1 
-1 + i 
1 







1 



1 + i 
1 



(21) 



and select an initial pure state of the system as \ip{0)) = |0) G T-L^. The state of the system at 
some specific time t is described by the transformed state \ip{t)) = U{t)\ip{0)) . 



e? 0.0 




» 0.0 




a 0.0 




(a) A, 



(b) 



(c) A, 



Figure 4. Shadow of the operators Ai, A2 and A3 of order three serves as a background 
for the trajectory representing the unitary dynamics defined hy U = exp{—iHt) with 



H given by (21 1 with the initial state 10(0)) (marked by circle in the picture). All 



plots are done for matrices translated in such a way that their trace is equal to zero 
and suitably rescaled. 



In order to use the numerical shadow to study the time evolution of the system, one 
needs to choose an arbitrary 3x3 non-hermitian matrix. To get some information on the 
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dynamics in the space of pure states of a qutrit and to observe it from different points of view, 
we selected the fohowing matrices 



A, 
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,^3 = 
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(22) 



For each of these auxihary matrices, the quantum dynamics can be now visuahzed as a 
trajectory in the complex plane defined by a parametric equation, 

z{t) := {m\A\m) = {m\A'\m), 

where the unitarily transformed matrix reads A' = U*AU 
of the initial state IV'(O)) = (1,0,0)^ G T-L^. generated by the Hamiltonian (21) is shown in 
Fig. |4]from three different perspectives determined by matrices (22). 



Ae 



-iHt 



(23) 

The time evolution 



5.1. Identical trajectories 

For given matrix A and hermitian matrix H a unitary time evolution induces a path in the 
numerical range given by 

Ce'^^^Ae'^^C, for given starting point ^. (24) 

In the mixed state scenario the trajectory is given by 

tr /3e~'^*^e^^*, for given starting point p. (25) 

The question one may pose is: under what conditions for two different starting points po ^ind 
pi trajectories on numerical range of A are identical: 

trpoe-'^'Ae'^' = tipie-'"'Ae'^\ teR. (26) 

To convince yourself that such a situation may occur consider the daily rotation of the earth 
around its axis. Choosing two initial points at the same meridian on opposite sides of the 
equator (say close to Cairo and Durban in Africa), we see that the trajectories they generate 
after projecting onto the equatorial plane do coincide. This is because the dynamics, both 
initial points and the kind of the projection, are chosen in a special way and satisfy certain 
constraints. To characterize these constrains in a general setting we start with the following 
definitions. 

For given matrix A, let Xa = {B £ M^(C) : trS = 0, ti BA = 0}. We also define the 
set Ha 

Ha = {H e (C) : Vt > 0, 5 G Xa we have M^iHt{B) G Xa], (27) 

where Ad is the adjoint mapping given by Kdc{B) = CBC^^. 
Now we can state fact concerning identical trajectories. 

Lemma 7 Trajectories tr pQe~^^^ Ae^^^ and tr pie~'^*yle'^* for t G M are identical if and only 
if Po — Pi ^ and H G Ha ■ 

The definition of Ha is somehow complicated, here we put reasoning which presents it 
in a simpler form. We have following property Adc = e^'^^ , where adc(i?) = [C, i?] (see 
e.g. [34J). Using above we can state the following lemma. 

Lemma 8 Hermitian matrix H is an element of Ha, if and only if, for all B G Xa we have 
adiH{B) G Xa. 
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Proof. If adiH{B) G Xa for all B £ Xa, then by iterating we have that ad^jj{B) G Xa for 
k = 0,1, . . .. Since Xa is a linear space we obtain, that 

OO 

M,MB) = ^H'i/l^) e Xa. (28) 

fc=0 

On the other hand if Ad^itH{B) G Xa, then using the fact that 

piHt TDp-iHt _ D 

■i^-^i = !To 1 • (^») 

and the continuity of the function X i— )■ trX^l we obtain the result. □ 

Note, that the condition i[H, B] G Xa can be stated as ti H[A, B] = 0, this follows from 
the cyclicity of trace. The linear space Xa is a real A^^ — 1 — d{A) dimensional space, where 
d{A) = dim({3f?(^),9(^)}), where ^{A) = \{A + A*) and 9(^) = ^(^ - A*). Thus in the 
generic case Xa has dimension A^^ — 3. The set Ha forms a real subspace of hermitian matrices 
orthogonal to sum of two real subspaces (adiSR(^)(XA) and adiQ(^)(X^))), 

Ha = (adij5(^)(X^) U ade(A)(^A))^, (30) 
where denotes the orthogonal component in the real space of hermitian matrices. 

6. Mixed states numerical shadow 

The standard numerical shadow ([2]) of matrix A is defined by choosing randomly a pure 
state p = \ip){ip\ with respect to the unitarily invariant, natural measure on the set of pure 
states, and taking the expectation value Tr^p. However, one may also consider an expression 
analogous to ^ and use it with a different measure in the set Qat of mixed states. More 
precisely, we introduce the mixed states numerical shadow of A with respect to a measure fi, 



P^{z):= / dp{p)6[z-TrpA . (31) 
JQn ^ ' 

The measure p. defined on the set Qat of mixed states of size A'" is supposed to be unitarily 
invariant. For instance, we will use the family of induced measures px obtained by taking 
a random pure state |.^) G Hn ® 'Hr and generating a mixed state by partial trace over the 
K-dimensional subsystem, p = Trx|0('?l- Since the pure states |^) are generated randomly 
the unitary matrices determining the eigenvectors of p are distributed according to the Haar 
measure on U{N). The probability distribution of the eigenvalues Aj of the random mixed 
state p of size N obtained in this way reads 

/ N \ N 

PN,K{^)=CN,Ksil-Y,^^] n^f"'^n(^^-^^)'- (32) 
\ i=l / i=\ i<j 

It is assumed here that K > N and the normalization constants Cn,k are given in [27]. In 
the symmetric case, K = N, the above formula simplifies and the measure pN coincides with 
the flat Hilbert-Schmidt measure, induced by the metric Q. In the opposite case K < N, 



the joint probability density function is given by (32) with exchanged parameters N ^ K. 

Consider now a pure state |.^) on the bi-partite N x K system. It can be represented in 
its Schmidt decomposition [2], 

mm{N,K} 



\0= E (33) 



i=l 
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where {|ej)}^^ is an orthonormal basis in T-L^ while is an orthonormal basis oH-Lk- 

Taking a partial trace of the projector over the K dimensional system we see that the 

spectrum of the resulting mixed state p coincides with the set of the Schmidt coefficients {Aj} 
of the pure state Thus formula (32 ) describes the distribution of the Schmidt coefficients of 
a pure state \S) drawn randomly according to the uniform distribution on the sphere S^^~^. 
By construction of the Schmidt decomposition of a random state |^), the vectors le,) and \ fi) 
can be considered as columns of unitary matrix in U{N) and U{K) respectively, distributed 
according to the Haar measure on the unitary group. A simple calculation shows that 

ra\n{N,K} 

min{7V,X} iaiii{N,K} 

= A,(e,|^|ei)= Ai(C/Uc/)i,i, (34) 

i=l i=l 

where [/ is a unitary matrix distributed according to the Haar measure on U{N). 

These considerations imply that the shadow oi A®^k is a mixture of diagonal elements 
of ^ in a random basis, given by the sum Y.T=i^'^^ ^iiU^^U\i. As before U stands for 
a random unitary matrix of size while the joint probability distribution function of the 
coefficients Aj is given by (32). This proves that the mixed states shadow of A with respect 
to the induced measure p,K coincides with the standard numerical shadow of the extended 
operator A® Ik-, 

PTi^) = PAmA^) = Pik®a{z). (35) 

The last equality follows from the property ( |20[ ). In the most important case, K = N, the 
induced measure fiN is equivalent to the Euclidean (flat) measure in corresponding to 

the Hilbert-Schmidt distance (|4]). Thus the projection of the 'full' set Qat of mixed quantum 
states on the plane determined by a given matrix A of size A^ is equivalent to the standard 
shadow of an extended operator ^ (8) Ijy. In the case of A^ = 2 this is visualized in figure [5|:), 
in which the shadow of the full Bloch ball Q2 can be compared with the shadow of the hollow 
Bloch sphere = 5^, displayed in figure [Hja). 

Note that for A' = 1 the induced measure fii is supported on the set 17 at of pure states 
only and coincides with the Fubini-Study measure, so formula (31) with // = /xi reduces to 
the standard definition ([2]) of the pure states numerical shadow. 



7. Large limit and random matrices 

It is instructive to analyse the numerical shadow of a random matrix in the limit of large 
matrix dimension A^. Let us consider two cases of the problem: the shadow of a random 
density matrix a generated according to the induced measures [27] and the shadow of random 
unitary matrix U distributed with respect to the Haar measure on U{N). 

The numerical shadow of random matrix which is distributed with unitarily invariant 
measure is related to the distribution of its arbitrary diagonal element in a fixed basis. In 
this section we consider the measures induced by partial trace and the Haar measure on the 
unitary group, which are unitarily invariant. Let us start with the following 

Lemma 9 Let A be a random square matrix of order N distributed according to a unitarily 
invariant measure. Let \x) be a random pure state of size N generated according to the Fubini- 
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Figure 5. Sketch of a projections onto a two-plane of a) ttie set C4 — A3 of — 4 
classical states onto the quadrangle formed by the numerical range W of a, normal 
matrix A of order 4; b) the set 5^2 = S'^ of one-qubit pure quantum states onto a 
disk formed by the numerical range of a non-normal Jordan matrix J2 of order 2; c) 
a mixed-states numerical shadow of J2, corresponding to the projection of the full 
3-dimensional Bloch ball onto a plane, is equal to the standard, pure states shadow of 
an extended matrix J2 18) l2- The picture is plotted using perspective. 



Study measure otlVLn = CP^ ^. Then the expectation value has the same distribution as the 
matrix element Ai^i 

P {{x\A\x)) = P {Ai^i) . (36) 

Proof. Since |x) is a random pure state, thus \x) ~ f^|0), where |0) is an arbitrary fixed state 
while f/ is a random unitary matrix of size A^. Now we write 

P ((x|^|x)) = P (^(0|[/U?7|0)) = P . (37) 

The last equality follows from invariance of the distribution P{A) with respect to unitary 
transformations. □ 

7.1. Shadow of random quantum state 

Consider a random density matrix a of size N generated with respect to the Hilbert-Schmidt 



measure ^hSi so that its eigenvalues Aj are distributed according to eq. (32) with K = N [27] . 
The diagonal elements of a are of the form 



where ^ij and r/,j are independent, identically distributed random variables with normal 
distribution A/'(0, 1). Basic properties of the Gamma distribution r(a, b) |28] imply that 

G 
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where Gi and G2 are stochastically independent variables distributed according to the Gamma 
distribution r(A^, 2) and T{N{N — 1),2) respectively. Therefore the diagonal elements of 
a random matrix a generated according to the measure ^hs are described by the Beta 
distribution with parameters {N,N{N — 1)}. 

The same reasoning can also be used for a general class of induced measures (32) 
parameterized by the size K of the auxiliary subsystem. In this case the diagonal elements of a 
density matrix a £ Qn generated with respect to the measure fj.N,K are distributed according 
to the Beta distribution with parameters {K,K{N — 1)}. 

Using the above reasoning and Lemma [9] we get the following 

Lemma 10 The numerical shadow of a random matrix a generated with respect to the induced 
measure ^n.k is given by the Beta distribution with parameters {K, K{N — 1)}, which can be 
expressed in terms of the Beta function, 



7.2. Random unitary matrices 

Let us now consider a random unitary matrix U distributed according to the Haar measure. 

Lemma 11 The numerical shadow of a Haar random unitary matrix is supported in the unit 
disk. This distribution is invariant with respect to rotations, and \{x\U\x)\'^ is distributed 
according to the Beta distribution with parameters {1, — 1}. 

Proof. Random unitary matrix distributed with Haar measure can be generated using the QR 
decomposition of matrices pertaining to the Ginibre ensemble [29] . The QR factorization can 
be realized by a Gram-Schmidt orthogonalization procedure. Then the element Ui^i of the 
generated unitary matrix reads 

= , ^^'^ , (41) 



where ^ is a nonhermitian random matrix from the Ginibre ensemble. Therefore 

Z]i=i(€i\i + 'nil) 



where and r]ij are independent, identically distributed random variables with normal 
distribution A/'(0, 1). Thus |C/i,iP has the Beta distribution with parameters {1, — 1} and 
using Lemma [9] we arrive at the desired result. □ 



8. Concluding remarks 

Our study may be briefly summarized by the following observation. The numerical shadow of 
a normal operator acting on T^jv reflects the structure of the set of (mixed) classical states, 
which belong to the probability simplex A^r.i, while investigation of numerical shadows of 
non-normal operators provides information about the set Qn of quantum states of size A^. 

In particular we have shown that the set of orthogonal projections of the set Qat of 
density matrices onto a two-plane is equivalent, up to shift and rescaling, to the set of all 
possible numerical ranges ll^(^) of matrices of order N. The numerical shadow of A forms a 
probability distribution on the plane, supported in VF(^), which corresponds to the 'shadow' 
of the complex projective space CP^~^ covered uniformly according to the Fubini-Study 
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measure, and projected onto the plane. Another probabiHty distribution in is obtained 

if one projects onto this plane entire convex set Qn of density matrices. If this set is covered 
uniformly with respect to the Hilbert-Schmidt (Euclidean) measure, an explicit expression for 
this distribution is derived. In this way the analysis of numerical ranges and numerical shadows 
of matrices of a fixed size N contributes to our understanding of the intricate geometry of the 
set Qat of quantum states [2j. 

The numerical range |3^0J and its generalizations [STJ [TC] found several application in 
various problems of quantum information theory. In analogy to the product numerical range, 
defined for spaces with a tensor product structure |32] , one can introduce the numerical shadow 
restricted to the subset of separable (product) states or the set of maximally entangled states 
|12j . Analysing such restricted numerical shadows for operators of a composite dimension NM 
one may thus investigate the geometry of the selected set of separable (maximally entangled) 
quantum pure states. Such an approach is advocated in a forthcoming publication [33] . 
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